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Abstract 



We present an evaluation of the non-mesonic decay widths for A-hypernuclei 
{AN NN, ANN NNN) within the framework of the polarization 
propagator method. The full A self-energy is evaluated microscopically in 
nuclear matter by using the functional approach, which supplies a theoreti- 
cally well grounded approximation scheme for the classification of the relevant 
diagrams, according to the prescriptions of the bosonic loop expansion. We 
employ average Fermi momenta, suitably adapted to different mass number 
regions (medium- light, medium and heavy hypernuclei) . Moreover, we study 
the dependence of the decay rates on the A^A^ and AA^ short range correla- 
tions. With a proper choice of the parameters which control these correlations 
in the new approximation scheme, it is possible to reproduce the experimental 
decay widths for yl > 10 hypernuclei. 

21.80.-ha, 13.75.Ev, 25.40.-h 
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I. INTRODUCTION 



Among the existing hypernuclei, those which contain one A hyperon are the most stable 
with respect to the strong interaction and they are the subject of this paper. The study 
of hypernuclear physics is connected to various and general problems of both nuclear and 
particle physics. In fact, it helps in understanding important aspects of weak interactions in 
nuclei, the role of sub-nucleonic degrees of freedom in the hadron-hadron interactions and 
the renormalization properties of mesons and hyperons in the nuclear medium. 

The most interesting hypernuclear decays are those involving weak processes, which 
directly concern the hyperon. The subject of the weak decay rates of hyperons embedded 
in nuclei has a quite long history, both on the theoretical and the experimental sides, but it 
has received a broader attention only in the last ten years. The weak decay of A-hypernuclei 
occurs through two different modes: the so called mesonic channel: 

A^ttN {Tm), (1.1) 

and the non-mesonic one, which can be mainly attributed to the following processes: 

AA^^A^A^ (Fi), (1.2) 

AA^A^^A^A^A^ (Fa), (1.3) 

both mediated by the exchange of n, p,rj,uj, K, K*, etc. The non-mesonic mode is only 
possible in nuclear systems and, nowadays, due to the difficulty of using A beams, the study 
of the A decay in nuclei is the only way to get information on the weak process AA^ NN. 

The free A decay is compatible with the A J = 1/2 isospin rule, which is based on the 
observed ratio FA^Tr-p/rA^TrOn — 2, but it is not yet understood on theoretical grounds. 
Besides, because of nuclear shell effects, it is not yet clear the level of violation of this rule 
for the mesonic decay in nuclear systems. In the present calculation we will assume this rule 
as valid. The momentum of the outgoing nucleon in the process A ttN is about 100 MeV, 
hence this decay mode is suppressed by the Pauli principle in nuclei, particularly in heavy 
systems (in infinite nuclear matter being strictly forbidden): the mesonic width is found to 
decrease fast as the mass number A of the hypernucleus increases This is confirmed 

by the few available experimental data. 

In the non-mesonic process AA^ — A^A^, the final nucleons have large momenta (~ 420 
MeV) and, apart from the s-shell hypernuclei, it dominates over the mesonic decay. The 
presence of large momentum transfers in the non-mesonic channels implies that the details 
of the nuclear structure have only little infiuence on the decay, but, on the other hand, 
the A^A^ and AA^ short range correlations are crucial. There appears to be an anticorre- 
lation between mesonic and non-mesonic decay, such that the total width is quite stable 
from light to heavy hypernuclei. In fact, as discussed in ref. the non-mesonic rates Fi 
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and r2 increase and reach saturation values as the mass number of the nucleus increases: 
r^M(rPb)/r^A/(f C) ~ 1.4, where T^m = Ti + T^. 

In this paper we present an evaluation of the decay rates for a A in nuclear matter within 
the framework of the polarization propagators. The relevant Feynman diagrams have been 
obtained following a functional approach, which allows to divide these diagrams into classes, 
according to the prescription of the so-called bosonic loop expansion (BLE). The calculation 
in finite nuclei (using the local density approximation) is not possible here because of the long 
computing time already for the calculation of the decay widths at fixed Fermi momentum 
(namely in nuclear matter). Nevertheless, in order to consider different mass regions, we 
assign an average Fermi momentum felt by the hyperon in the various nuclei. 

We remind that the baryon-baryon strong interactions cannot be treated with the stan- 
dard perturbative method. However, in the study of nuclear phenomena we always need to 
sum up the relevant diagrams. For instance,one usually performs the summation of the infi- 
nite classes of diagrams entailed by the RPA and Hartree-Fock approximations. However, in 
the above quoted schemes no prescription is given to evaluate the "next-to-leading" order. 

The functional techniques can provide a theoretically founded derivation of new classes 
of expansion in terms of powers of suitably chosen parameters. On the other hand, as we will 
see, the RPA automatically appears in this framework as the mean field level. The method 
has been extensively applied to different processes in nuclear physics [^-0. We will use it 
for the calculation of the A decay rates in nuclear matter, which can be expressed through 
the nuclear responses to pseudoscalar-isovector and vector-isovector fields. We will see that 
the responses include ring-dressed meson propagators (which represent the mean field level 
of the theory) and almost the whole spectrum of 2p-2h excitations (expressed in terms of a 
one-loop expansion with respect to the ring-dressed meson propagators) which are required 
for the evaluation of Actually, the semiclassical expansion leads to the prescription 
of grouping the relevant Feynmam diagrams in a consistent many-body description of the 
"in medium" meson self-energies: the general theorems and sum rules of the theory are 
preserved. 

The paper is organized as follows. In Sec. |l| we summarize the model used for the 
calculation of the decay rates. In Sec. |lll| the functional approach to the spin-isospin nuclear 
response functions in presented. By using the semiclassical approximation we will introduce 
the class of Feynman diagrams needed for the evaluation of Our results are presented 
and discussed in Sec. |3 we first study the sensitivity of the decay rates to the NN and AA^ 
short range correlations, and then we parametrize these correlations in order to reproduce 
the experimental decay widths for three different mass regions in the hypernuclear spectrum. 
Our conclusions are given in Sec. 0. 
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II. POLARIZATION PROPAGATOR METHOD 



In this section we briefly summarize how the A decay in nuclear matter can be studied 
employing the polarization propagator method 0J|,||]. This technique provides a unified 
picture of the different decay channels and it is equivalent to the standard wave function 
method, used by other authors in refs. P,p|-pT|,p|. In a previous work |1| the evaluation of 
the decay widths was first performed in nuclear matter and then extended to finite nuclei 
via the local density approximation: in this framework the experimental partial decay rates 
in a range of nuclei from to j^^Ph have been reproduced. The two-body induced decay 
ANN —>■ NNN has been accounted for via a purely phenomenological parametrization 
(using data on pionic atoms) of the 2p-2h polarization propagator which enters the A self- 
energy Sa. The total decay rate is obtained through the relation: 

Ta = -2ImSA. (2.1) 

Here we present a microscopic calculation of the 2p-2h polarization propagator within the 
theoretically consistent scheme which will be introduced in the next section. 

Let us first remind the main steps in the evaluation of Fa. We start from the A — > vrA^ 
effective lagrangian: 

Ca^n = Gml'^j^iA + B-i^)T ■ cj^^^A + h.c, (2.2) 

where the values of the weak coupling constants G ~ 2.211 ■ 10^' /m"^, A = 1.06, B = —7.10 
are fixed from the free A decay. The constants A and B determine the strengths of the parity 
violating and parity conserving A — > vrA^ amplitudes, respectively. In order to enforce the 
A/ = 1/2 rule, in eq. ( |2.2|) the hyperon is assumed to be an isospin spurion with = —1/2. 
Using the Feynman rules, in the non relativistic limit the A self-energy in nuclear matter is 
(see fig. |1|): 

SA(fc) = 3^{Gmlf | ^ js^ + ^g^j ^.^(^^^^(^ _ (2.3) 

Here, S = A, P = rriT^B /2mN, while: 

0{\p\-i 

Po - En{p) -Vn + ie ' po- Ei^{p) -Vn -ie 

and: 



G.M^ "'1^1-^'. + fr'"'' (2.4) 



GT,{q) — — ^ ^ / (2-5) 

are the nucleon and pion propagators in nuclear matter, respectively. In the above, p = 
{po,p) and q = {qo,q) denote four-vectors, kp is the Fermi momentum, is the nucleon 
total free energy, Vn is the nucleon binding energy (assumed constant), and S* is the pion 
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FIG. 1. A self energy in nuclear matter (the dashed line represents a dressed pion). 



proper self-energy. Moreover, in eq. ( |2.3|) we have included a monopole form factor for the 
vrAA^ vertex: 

FAq) = ^ (2.6) 

with the same cut-off = 1.3 GeV used for the tcNN strong vertex. In fig. |^ we explicitely 
show the lowest order Feynman graphs for the A self-energy. Diagram (a) represents the 
bare self-energy term, including the effects of Pauli principle and of the binding on the 
intermediate nucleon. In (b) and (c) the pion couples to a particle-hole (p-h) and a A-h 
pairs, respectively. In diagrams (d) and (e) we show examples of 2p-2h excitations, while (f ) 
is a RPA iteration of diagram (b). Once evaluated the integral over go in ( p.3|) , the nuclear 
matter A decay width (eq. ( p.l|) ) becomes [^: 

TA{k, p) = -6{Gmiy J j^Oilk -q\- kF)e{k, - E^{k - q) - V^) 

xlma(g) \q^^=ko-Et,(k-q)-VM^ (2-7) 

where 



l-VLiqpLiq) l-VTiqWriqV 

consist of a longitudinal and a transverse part. In eq. ( p2.7| ) the first 6'-function forbids 
intermediate nucleon momenta (see fig. |1|) smaller than the Fermi momentum and the second 
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(d) (e) (f) 

FIG. 2. Lowest order terms for the A self-energy in nuclear matter. The meaning of the various 
diagramms is explained in the text. 
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one requires the pion energy go to be positive. Moreover, the A energy, ko = E^^k) + V\, 
contains a binding term. The pion hnes of fig. |^ have been replaced in eq. ( p.8| ) by the 
interactions Vl, Vt, S, Pl, Pt, which include vr and p exchange modulated by the effect of 
short range correlations. They have the following expressions: 



V,{q) = ^{q'FMGUq)+9L{q)} 



(2.9) 



VT{q) = ^ {q'C,Fl{q)Gl{q)+gT{q)] . 



(2.10) 



{q'F'^{q)Gl{q)+gt{q)}, 



(2.11) 



pTiq) = ——9^{q), 



(2.12) 



~S{q) = {F^{q)Gl{q) - F^{q)Gliq)} \ q 



(2.13) 



In the above, G%j = \/{ql — q^ — m\j) denotes a free meson propagator and Fp is the pNN 
form factor (eq. ( |2.6| ) with cut-off Ap = 2.5 GeV). The exchange of p-mesons is brought 
into play by the short range correlation Pt embodied in the AA^ NN potential. Form 
factors and propagators with a tilde imply that they are calculated with the replacement 
qr^ — > q^ + q^i where 1/gc is a lenght of the order of the hard core radius of the interaction. 



and Gp is the ratio: 



a, 



P 

rni 



J TT 



(2.14) 



The short range correlations included in the potentials of eqs. ( |2.9| )- (|2T^) are explicitely 
reported in the appendix of ref . . Their momentum dependence (governed by the param- 
eter gc) ensures the convergence of the integrals contained in diagrams like (d) and (e) of 
fig. (this behaviour is not crucial in ring approximation). In the calculation, the values of 
the correlation functions in g = 0, 



9' = 9lA^)^ 9 a 



9tA0) 



(2.15) 



(namely the Landau parameters), will be considered as free. 

Furthermore, in eq. ( p. 81) Ul and Ut contain the Lindhard functions for p-h and A-h 



excitations |jT2| and also account for the irriducible 2p-2h polarization propagator: 

UL,T{q) = u'\q) + u^\q) + ul^?{q)- 



(2.16) 
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They appear in eq. (|2.8|) within the standard RPA expression. The Lindhard functions U^^, 
U^^ are normahzed as in ref . . In the following section we will develop an approximation 



scheme for the classification of the various diagrams to be included in the calculation of the 
A self-energy. Eq. ( |2.7| ) depends, both explicitly and through Ul,t{(i), on the nuclear matter 
density p = 2fc|./37r^. 

We remind that in nuclear matter the mesonic decay is strictly forbidden, the final 
nucleon in A ^ ttN having a momentum lower than the Fermi one (for a A decaying at 
rest, q ^ 100 MeV < kp c::^ 270 MeV). In finite nuclei the mesonic rate is dominant only in 
s-shell hypernuclei and rapidly decreases with the mass number. 

III. FUNCTIONAL APPROACH TO THE SPIN-ISOSPIN NUCLEAR 

RESPONSE FUNCTIONS 

We first evaluate the polarization propagator in the pionic (spin-longitudinal) channel. 
In order to exemplify we consider a lagrangian describing a system of nucleons interacting 
with pions through a pseudoscalar coupling: 

C^r^ = ^{i^ - m^)^ + ]^d^(j) ■ d^(j) - ^mlcjy^ - i^T^ ■ cji, (3.1) 

where if) {(p) is the nucleonic (pionic) field, and: 

r = ^75T (3.2) 

{g = 2/7rmAr/m,r) is the spin-isospin matrix in the spin-longitudinal isovector channel. We 
remind that in the calculation of the hypernuclear decay rates we also need the polarization 
propagator in the transverse channel (see eqs. ( |2.7| ),( p^ ): hence, we will have to include in 
the model another mesonic degree of freedom, the p meson. As we will see, the semiclassical 
expansion is characterized by the topology of the diagrams, so the same scheme can be 
applied to mesonic fields other than the pionic one. In this section we use a relativistic 
formalism, the non-relativistic reduction of our results being trivial. 

Let us now introduce a classical external field with the quantum numbers of the pion. 
The lagrangian then becomes: 

C^N C^N - ii'Ttp ■ (p. (3.3) 
The corresponding generating functional in terms of Feynman path integrals has the form: 

Z[(p] = j V ip,ip,(f) exp |i J dx Ct,n{x) — iip{x)Tip{x) ■ (p{x) | (3.4) 

(here and in the following the coordinate integrals are 4-dimensional). All the fields in 
the functional integrals have to be considered as classical variables, but with the correct 
commuting properties. The physical quantities of interest for the problem are deduced from 
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the generating functional by means of functional differentiations. By introducing a new 
functional such that 



Z\^\ = exp{iZc[ifi]}, 



(3.5) 



the polarization propagator turn out to be the second functional derivative of Zc with respect 
to the source ip of the pionic field: 



Uij{x,y) = - 



(3.6) 



We notice that the use of Z^. instead of Z in eq. ( |3.6[ ) amounts to cancel the disconnected 
diagrams of the perturbative expansion (linked cluster theorem). From the generating func- 
tional Z one can obtain different approximation schemes according to the order in which 
the functional integrations are performed. By integrating eq. ( p.4| ) over the mesonic degrees 
of freedom first, the generating functional can be written in terms of a fermionic effective 
action S^q. Up to an irrelevant multiplicative constant: 



Z[p] 



V 



(3.7) 



The remaining integration variables are interpreted as physical fields and S*^ describes a 
quadrilinear non-local time- or energy- dependent nucleon-nucleon interaction induced by the 
exchange of one pion: 



= / dxdy [ij{x) G,^^{x - y)ij{y) 



+- J2 ij{x)T,ij{x)Gl{x - y) ^{y)T,ij{y) 
^ i=i 



(3.8) 



where Gn and G° are the nucleon and free pion propagators, respectively, which satisfy the 
following field equations: 



rriN -iT ■ (p) Gn{x - y) = S{x - y), 



(3.9) 



+ ml) Glix -y) = -5{x - y). (3.10) 

The pion propagator is diagonal in the isospin indices: (G°)^^. = 5ijG^. The effective action 
( p.8|) can then be utilized in the framework of ordinary perturbation theory and will not be 
employed in the following. In fact, it does not bring significant novelties with respect to the 
usual calculations; furthermore, it cannot be correctly renormalized due to the absence of a 
term proportional to (/)^, which is needed to cancel the divergence of the 4-points fermion 
loops. 
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A. The bosonic effective action 



Alternatively, it is possible to eliminate in eq. (|3.4| ) the nucleonic degrees of freedom 
first (without destroying the renormalizability of the theory [^]). Introducing the change of 
variable (f) ^ <p — cp, eq. (|3.4| ) becomes: 



Z[ip] = expj^y dxdy(p{x) ■ \x - y)(fi{y) 



X / V 



ip,ip,(f) exp <^ i / dxdy ■il){x)G {x — y)il){y) 



+ -ct>[x)-Gl-\x-y) (0(y) + 2(^(2/)) 



where the integral over %p,%p 



IS gaussian: 



V 



i/j^ip expli / dxdyil){x)Gj^{x~y)il){y)\ = (detCjAr 



(3.11) 



(3.12) 



Hence, multiplying eq. (|3.11| ) by the unessential factor det (G^ being the free nucleon 
propagator), which only redefines the normalization constant of the generating functional, 
and using the property detX = exp {Tr InX}, we obtain: 

Z[ip] = exp 1^ I dxdy v{x) ■ Gl\x - y^y)^ J V [0] exp [iS^^ [0]} , (3.13) 



with: 



^fff [0] 



dxdy \ ^-ct>{x) ■ Gl\x - y) [ct>{y) + 2^{y)] + K[0] \ , 



(3.14) 



OO 1 



n=l 



n 



i^Tr(rirj) j dx dy Yi^{x,y)(j)i{x)(j)j{y) 

+^ E Tr (r^r.-r,,) j dxdydzU%x,y,z)^i{x)^j{y)Mz) + 0{<l)' 



where 



-tlf{x,y) = tGUx-y)tG%iy-x), 
-in°(x, y, z) = iG%{x - y)iG%{y - z)iG%{z - x), 



etc. 



(3.15) 



(3.16) 
(3.17) 



We have thus derived an effective action for the bosonic field </>. This action contains a term 
for the free pion field and also a highly non-local pion self-interaction VJr, which is described 
by the Feynman diagrams shown in fig. ^ This effective interaction is given by the sum of 
all the diagrams containing one closed fermion loop. We note that the function of eq. ( p.l6| ) 
is the free polarization propagator, namely the Lindhard function. Moreover, the functions 
Il'^{x, y, ..,z) are symmetric for cyclic permutations of the arguments. 
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FIG. 3. Diagrammatic representation of the bosonic effective action ( |3.14 ) 



B. Semiclassical expansion 

The next step is the evaluation of the functional integral over the bosonic degrees of 
freedom in eq. ( |3.13|) . A perturbative approach to the bosonic effective action (|3.14|) does 



not seem to provide any valuable results within the capabilities of the present computing 
tools and we will use here another approximation scheme, namely the semiclassical method. 

Mean field level 

The lowest order of the semiclassical expansion is the stationary phase approximation 
(also called saddle point approximation (SPA) in the euclidean space): we require the bosonic 
effective action to be stationary with respect to arbitrary variations of the fields (pf 

d(pi{x) 



From the partial derivative of eq. (|3.14 ) we obtain for the classical field (f) the following 
equation of motion: 



(□ + m^)0,(x)= fdyG:-\x~y)^,{y) + -£fl, (3.19) 



SV. [0] 
(50i(x) 

whose solutions are functional of the external source (f. The exact solution cannot be written 
down explicitly. However, due to the particular form of when (p ^ 0, one solution is 

= 0; we can then express the general solution of eq. (|3.19| ) as an expansion in powers of 

= E / dyAj{x,y)ipjiy) + / dy dz Bijk{x,y, z)ipj{y)Mz) + O (y^) . (3.20) 

3 ^ j,k •' 



Substituting eqs. ( p.20| ) and ( |3.15| ) into ( p.l9| ) and keeping the terms linear in yjj we obtain 
the following relation for Aif 

Aij{x,y) - Tr (f^) j dudvGl{x - u)U%u,v)Aij{v,y) = 5ij5{x - y). (3.21) 

Finally, by introducing the ring-dressed pion propagator G^^^^, which satisfies the Dyson 
equation: 
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^(x -y) = Gl{x - y) + Tr (r^) j dudvGl{x - u)Ii\u, v)GT^{v - y), 
or, formally: 

the solution of (|3.21| ) reads: 



1 - Tr (rf) Gono 



(3.22) 



(3.23) 



(3.24) 



A,,(x,y) = 6,, j dzGT^{x - z)GV{z - y). 
Therefore, the saddle point of ( ^.141 ) at first order in the source ip is: 

0r(x) = j dydzGT^{x - z)Gl~\z - y)^,{y) ^ J dy (Cf ^G^') {x - y)^,{y), (3.25) 
and the corresponding bosonic effective action reads: 







ring 



- / dxdydudvGl' {x - u)^{u) ■ G'^^'^ix - y)Gl~ {y-v)^{v). (3.26) 



Now, the generating functional of eq. ( p.l3| ) takes the form: 



Z = exp <j ^ I dxdydu dv ^{u) ■ G° \x — u) 



X 



G^ix - y) - GT'ix - y)] Gl-\y - v^v)] 



(3.27) 



The polarization propagator is then evaluated using eqs. (|3.5|), (|3.6|). We obtain that in the 
saddle point approximation it coincides with the well known ring expression: 



^ij{x,y) = 5, 



n°(x, y) + Tr (rn / du dv U%x, M)Gfs(M - v)lf{v, y) 



or, formally: 



n 



no 



n 



ring 



(3.28) 



(3.29) 



1 - Tr (r?) GIW 

Hence, the ring approximation corresponds to the mean field level of the present effective 
theory. 



Quantum fluctuations around the mean field solution (one-loop corrections) 

In the next step of our semiclassical expansion we write the bosonic effective action as: 



S!^[4>\ = S!^[4P]+-Y.j dxdy 



S'S^^ [0] 
S(l>i{x)S(l)j{y)\^^^o 



0,(x)-0°(x)l U(y)-0O(y) 



(3.30) 
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where now (fp also contains the second order term in the source (see eq (|3.2CI|) ). Then, 
after performing the gaussian integration over 0, the generating functional reads: 



Z [ip] = exp <{ ^ / dxdycfiix) ■ G° \x - y)(p{y) 



X exp ^tSg [0° 

and the polarization propagator is: 

6^ 



— -Trln 



S'SS, [cf>] 
6(f)i{x)6(f)j{y) 



(3.31) 



cf>=4>o ^ 



6ipi{x)6ipj{y) 



'-'cS 



0° 



H — Tr In 
2 



^5(t)k{x)5(t)i{y) ^ 



(3.32) 



In the above the second derivative of the effective action (|3.14|) at the order 0^ turns out to 
be: 



5,,Gr (x-y) + Tr(r,r,)nO(x,y) 



6(f)i{x)6(f)j{y) 

+ ^ du \Tt [TiTjTk) n°(x, y, u) + Tr {TjTiTk) n°(y, x, u)] (f)k{u) 



+ / dudv \Tt {TiTjTkTi)Il^{x,y,u,v) + Tt {TjTiTkTi)Il°{y,x,u,v) 

k,l 

+ Tr (r.r^r.Tfc) Ii\x, v, y, u)] Mu)Uv)- (3.33) 



The second term in the r.h.s. of eq. ( |3.33| ) does not affect the calculation of eq. ( |3.32| ). By 
substituting ( p.20| ) in the equation of motion ( |3.19| ), from the terms of order we get the 
following Bijk functions: 



B,,,{x, y, z) = 2Tr (r.r.T,) / du dv dt U%u, v, t) 

xGT'ix - u) (GT'GV) {v - y) {GT'GV) it - z). 



(3.34) 



Now we have to calculate the logarithm in eq. ( p.32| ) up to second order in We can multiply 
the generating functional ( p.31[ ) by the factor (detG' 



inessential in the calculation of 



the polarization propagator (this corresponds to multiply eq. ( p.33|) by G°). Then, after 
calculating eq. ( |3.33| ) for cj) = 0°, with 0° given by the eqs. ( |3.2CI| ), ( p.24| ), ( p.34| ), we expand 
the logarithm up to (p^ and take the trace to the same order. Finally, the derivation with 
respect to the external source provides the following total polarization propagator: 



Iiij{x,y) = 6ijU{x,y), 



(3.35) 



with: 
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n(x, y) = n"'^s(x, 2/) + E Tr (r^r,) / du dv s(n - v)Il''{x, u, y, v) 

kl •' 

+ ^Tr(r,rO / dudvG^^^{u-v) [U''{x,u,v,y) + U''{x,y,v,u) 

kl 

+ J dudvdwdsG':^''^{u-w)G'^''^{v-s)U\x,u,v) 

X [Tr {TkTiTmTn) n'{y, w, s) + Tr (r,.r;r„r„) n°(y, s, w)] . (3.36) 



klmn 



We remind that the second derivative, with respect to the external source, of S^^ cf)"^^ 
and S^^ with 0"°^ [0°] given by eq. dOSj) [eqs. (lOOl) , ( F^ , (Hill)], gives the same 



result (the ring polarization propagator) when evaluated ad 95 = 0. 

The Feynman diagrams corresponding to eq. (|3.36|) are depicted in fig. ^. Diagram (a) 
represents the Lindhard function I\^{x,y), which is the first term of n"'^^(x, |/). In (b) we 
have an exchange diagram (the dashed lines represent ring-dressed pion propagators); (c) 
and (d) are self-energy diagrams, while in (e) and (f ) we show the correlation diagrams of our 
approximation. The approximation scheme we have developed is also referred to as bosonic 
loop expansion (BLE). The practical rule to classify the Feynman diagrams according to 
their order in the BLE is to reduce to a point all its fermionic lines and to count the number 
of bosonic loops left out. In our case the diagrams (b)-(f ) of fig. ^ reduce to a one-boson-loop. 
Diagrams (b), (c), (d) can be represented by the loop (A) of fig. ||, while (e), (f) correspond 
to the loop (B). 

When we include in the model the excitation of barionic resonances, then we have to 
replace the fermionic field with multiplets. The topology of the diagrams remains the same 
as in fig.^ but, introducing for example the A resonance, each fermionic line stands for a 
nucleon or for a A, taking care of isospin conservation. One thus obtain 15 exchange, 30 
self-energy and 98 correlation diagrams (see ref. ^ for the whole diagrammology) . 

Moreover, since this expansion is characterized by the topology of the diagrams, we can 
include in the model additional mesonic degrees of freedom together with phenomenological 
short range correlations. In particular, the extension to other spin-isospin channels simply 
amounts to change the definition of the Fj and the same occurs for the non-relativistic 
reduction of the theory. Accordingly, for the non-relativistic pion exchange Fj becomes 
(cr ■ q)Ti, for the p exchange it reads (cr x q)kTi, k being a spatial index, and for the u) 
exchange Fj = (cr x q)j. The exchange of w-mesons is taken into account only inside the 
one-boson-loop diagrams (b)-(f) of fig. ^, but not in the mesonic lines stemming from the A 
decay vertex, where the exchanged meson has to be an isovector (vr and p). Beyond vr, p and 
uj mesons, the present approach also contains (partly) the exchange of the scalar-isoscalar 
cr meson: as in the language of the Bonn potential [|l^], the latter is described through box 
diagrams (contained in the correlation diagrams of fig. namely by the exchange of two 
pions with the simultaneous excitation of one or two of the intermediate nucleons to a A. 

A further difficulty arises if we start from a potential theory, instead of one containing 
bosons as true degrees of freedom. This disease is however easily overcome by means of a 
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(e) 



(f) 



FIG. 4. Feynman diagrams for the polarization propagator of eq. ( ^.36 ): (a) particle-hole; (b) 
exchange; (c) and (d) self-energy-type; (e) and (f) correlation diagrams. Only the first contribution 
to the p-h ring expansion has been drawn. The dashed line represent ring-dressed pion propagators. 




(A) (B) 

FIG. 5. First order diagrams in the bosonic loop expansion. Diagrams (b), (c), (d) of fig. ^ 
reduce to diagram (A), while (e), (f) to (B). 
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Hubbard-Stratonovitch transformation, which enables us to substitute a potential with an 
interaction of the nucleons with a suitably introduced auxiliary field. For a scalar-isoscalar 
potential V the relevant identity reads: 

exp 1^ y" dxdyij{x)ij{x)V{x - y)i'{y)i^{y)^ = 

V det V J V[a]expl^-J dx dy a{x)V^^{x — y)o-{y) + i J dxilj{x)ip{x)a{x)'^ , (3.37) 

where a is the auxiliary field. Clearly, the previous derivation will remain valid, providing 
one substitutes the inverse propagator of the auxiliary field with the inverse potential in the 
"free" part of the action. 

Finally, a relevant point for the feasibility of the calculations is that all fermion loops in 
fig. ^ can be evaluated analytically [|l^,|l^], so that each diagram reduces to a 3-dimensional 
(numerical) integral. 

In particular, the formalism can be applied to evaluate both the longitudinal and the 
transverse functions Ul,t of eq. ( |2.16| ), which are required in eqs. (pT7|), (p.8|). In the 
approximation of fig. ^ we have to replace eq. ( |2.8| ) with: 



mr J ' ' l-Vi(g)f/i(g) 



where 



l-^L(g)f/i(g) l-VriqW^iq) 
+ IS'iq) + P|(g)l U^Aq) + ^l(g)f/2,T(g), (3.38) 



Ui = UP'' + [/^^ (3.39) 



and U2,L, f^2,T are evaluated from diagrams (b)-(f) of fig. § (taking into account the A 
excitation as well) using the standard Feynam rules. The normalization of these functions 
is such that U^^i^x.y) = 4n°(x,2/), with 11° given by eq. ( p.l6|) . 



IV. RESULTS AND DISCUSSION 

In the present work the evaluation of the hypernuclear decays is performed in nuclear 
matter. However, in order to compare with experimental data on finite (real) nuclei, we 
employ different Fermi momenta, which are fixed on the following basis. In local density 
approximation the "local" Fermi momentum kp{r) is related to the nuclear density: 



[a 



0.54 fm, R{A) = 1.12^1/3 - 0.86^-^3 f^)^ by the equation: 
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TABLE I. Average Fermi momenta. 







-pexp 






<kF> eq. Q 
(fm-i) 


< kp >eq. (4.4) 

(fm-i) 


Medium-Light: a B - f C 
Medium : f Si - fFe 
Heavy: f^Bi - f^U 


0.94 4 
1.20 
1.45 


- 1.07 
^ 1.30 
1.70 [ 


lie 
20- 


1 

22] 


0.95 
~ 1.1 
1.21 


1.08 
~ 1.2 
1.36 



Q -I 1/3 

2 / \ 

2^ 



(4.2) 



For the present purpose, an average, fixed Fermi momentum can be obtained by weighting 
each local kp with the density itself: 

1 f 

< kp >A= J J drkj{r)pA{r). (4.3) 

Alternatively (and more realistically) the average Fermi momentum should be determined 
by the probability density of the A inside the nucleus, according to the following definition: 

<kF>A= Jdrk^{r)\Mr)\', (4.4) 

ipA being the A wave function. We have calculated the latter from a Wood-Saxon potential 
with thickness a = 0.6 fm and radius and depth which reproduce the measured s and p 
A-levels ^ . Since the A wave function is preferably located in the interior of the nucleus (in 
fact, the hyperon occupies the Is level), we expect larger < kp > values from the prescription 



O than from (gj). 

We classify the hypernuclei for which are available experimental data on the non-mesonic 
decay rate into three mass regions (medium-light: A ^ 10; medium: A ^ 304-60; and heavy 
hypernuclei: A ~ 200), as shown in table |. The experimental bands include values of the 
non-mesonic widths which are compatible with all the quoted experiments. In the third and 
fourth columns we report the corresponding < kp >, as calculated with eqs. ( [4.3|) and ( [4. 41) , 
respectively. Accordingly, in the calculations of the decay width we use kp = 1.1 fm~^ for 
medium-light, kp = 1.2 fm~^ for medium and kp = 1.36 fm~^ for heavy hypernuclei. 

Before discussing the results, we list here the parameters which enter the calculation of 
the hypernuclear width. We point out that their values (with the exception of the Lan- 
dau parameters g' and g'/^, which will be discussed separately) have not been taken as free 
parameters. Instead, we kept them fixed on the basis of the existing phenomenology (for 
example in the analysis of the quasi-elastic electron-nucleus scattering, of the spin-isospin 
nuclear response functions, etc). 

• The correlation momentum is set to gc = 0.78 GeV; 
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The cut-offs in the pion, p and u form factors are: Attatat = 1.3 GeV, ^.■nNA = A,rAA = 
1 GeV, ApNN = ApATA = ApAA = 2.5 GeV, A^nn = Kaa = 1 GeV; 

The pion couphng constants are: /^^^/47r = 0.08, /^ata/^tt = 0.32, /^aa/^tt = 0.016; 

The rho couphng constants are (see eq. ( P^.14| )): CpAfA^ = CpAfA = CpAA = 2.3; 

The uj couphng constants are: C^nn = C<^aa = 1.5 (renormahzed values); 



• Finahy, the difference between the nucleon and A binding energies, which enter eq. 
as AV^ = Va — IV, has been fixed, for the heaviest nuclei {kp = 1.36 fm~^), from 
the corresponding depths of the binding potential (V/v = —55 MeV, Va = —32 MeV). 
Then, by assuming that AV is roughly proportional to the density, we obtain AV = 18 
MeV at kp = 1.2 fm and AV = 12 MeV at kp = 1.1 fm'^ It is worth pointing 
out, however, that the resulting widths are only weakly affected by the precise value 
of AV (a change of 5 MeV in AV would result in a 3% variation in the decay rates). 

The numerical evaluation of the first order contributions in the BLE [diagrams (b)-(f) 
in fig. Q presents a serious problem when the mesons (tt, p and uj in our case) are dynamic. 
The pionic branch in RPA is coupled to the particle-hole mode. By using non-relativistic 
kinematics for the nucleons, the relativistic pionic branch enters the particle-hole continuum 
for a momentum of about 1.6 GeV. This fact generates large oscillations in the calculation 
of the one- loop diagrams (where the momentum of the mesons is integrated out), which 
are absolutely unphysical and unavoidable from the numerical point of view. To prevent 
this mishap, we are forced to use static mesons in the evaluation of the one-boson-loop 
contributions. In all the mesonic lines of fig. ^ which are "external" to the bosonic loop we 
can use, on the contrary, the correct dynamical prescription, the relevant momenta being 
limited to about 0.6 GeV. 

An important ingredient in the calculation of the A decay rates is the short range part 
of the NN and AA^ interactions: in fact, in the non-mesonic processes the exchanged mo- 
mentum is very large (about 400 MeV). We parametrize these short range correlation with 
the functions reported in the appendix of ref. The zero energy and momentum limits 
of these correlations, g' and g'j^ (which we call, for hystorical reasons. Landau parameters), 
are considered as free parameters. No experimental constraints are available on g'^, while 
in the framework of the ring approximation (namely by neglecting the 2p-2h states in the 



A self-energy), realistic values of g' lie in the range 0.6 ^ 0.7 [^]. However, in the present 
context g' correlates not only p-h pairs, but also p-h with 2p-2h states and besides, in some 
diagrams [for instance (d) and (e) of fig. 0, two consecutive g' are connected to the same 
fermionic line, introducing some double counting, namely a renormalization of g'. In the 
picture of figures |^ and ^ the A self-energy acquires an energy and momentum behaviour 
which cannot be explained and simulated on the basis of the simple ring approximation. 
Therefore, the physical meaning of the Landau parameters is different in our scheme with 



18 



TABLE II. g' values compatible with the experiments. 





1-loop 


ring 


kp = 1.1 fm-i 


> 0.75 


0.45 ^0.65 


kp = 1.2 fm-i 


0.75 ^0.90 


0.55 0.65 


kp = 1.36 fm-i 


0.70 ^ 1.00 


0.65 ^0.75 



respect to the customary phenomenology, hence in the present paper we will use g' as free 
parameters, to be fixed in order to reproduce the experimental hypernuclear decay rates. 

As discussed in ref. |]1[, for fixed g' the non-mesonic width (the total width in nuclear 
matter, where Tm = 0) has a minimum as a function of g'/^, which is almost indipendent 
of the value of g'. This characteristic does not depend on the set of diagrams taken into 
account in the calculation, but it is simply due to the interplay between the longitudinal 
and transverse parts of the p-wave AA^ — > A^A^ potential [Pl and Pp functions of eqs. 
])]. In the present calculation the minimum is obtained, again, for g'^ ~ 0.4. 



Fixing g'j^ = 0.4, in ring approximation we can reproduce the experimental decay rates 
by using g' values which are compatible with the existing literature. In figure ^ we show, 
as a function of g' (for g'j^ = 0.4), the calculated non-mesonic decay widths (in unit of 
the free A width) for the three mass regions of table |. The thick solid curves refer to 
the one-boson-loop approximation of eq. (|3.38|) and fig. |^, while the dot-dashed curves are 



obtained through an RPA iteration of both the particle-hole and the one-boson-loop dia- 



grams, namely by using eq. (2^). However, we remind that only the former approximation 
has a theoretically founded basis, in line with the scheme introduced in sec. |T|: indeed, 
the present RPA calculation has the tendency to overestimate, in the acceptable range of 
g' values, the experimental non-mesonic widths. The dashed lines represent the pure ring 
approximation. The calculation is compatible with the experimental bands for the g' val- 
ues reported in table ^ As we have already noticed, the intervals corresponding to the 
ring calculation are in agreement with the phenomenology of other processes, like the (e, e') 
quasi-elastic scattering. However, only the full calculation (column 1-loop) allows for a good 
description (with the same g' values) of the rates in the whole range of kp considered here. 
In fig. 1^ we see the dependence of the non-mesonic widths on the Fermi momentum. The 
solid lines correspond to the one-loop approximation, with g' = 0.7, 0.8, 0.9 from the top to 
the bottom, while the dashed lines refer to the ring approximation, with g' = 0.5,0.6,0.7 
(again from the top to the bottom). We can then conclude that for the one- loop calculation 
the best choice of the Landau parameters is the following: 

g' = 0.8, g'^ = 0.4. (4.5) 

This parametrization is the same which was employed in a different theoretical framework 
M to reproduce the experimental decay rates in the range - ^^^Pb. However, we must 
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2.0 



1.5 



1.0 



(Medium-Light) 



kF=l.l fm" 



l-loop eq. (3.38) 
RPA eq. (2.8) 
ring 



EXP 



0.5 0.6 0.7 g' 0.8 0.9 
(Medium) 



2.0 



^ 1.5 ^ 



1.0 




0.5 0.6 0.7 g' 0.8 0.9 1 

(Heavy) 



2.0 



1.5 



1.0 




0.5 



0.6 



0.7 g' 0.8 



0.9 



FIG. 6. Dependence of the non-mesonic width on the Landau parameter (7', for g'^^ = 0.4. 
The three plots correspond to the classification of table |. The thick solid curves refer to the 
one-boson- loop approximation of eq. ( p.38| ), the dot-dashed ones to the RPA calculation of eq. ( |2.8D 
and the dashed ones to the ring approximation. The experimental bands of table | lies in between 
the horizontal solid lines. 
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1.8 
1.6 
1.4 
1.2 
1.0 




FIG. 7. Dependence of the non-mesonic width on the Fermi momentum. The sohd curves refer 
to the one-boson-Ioop approximation (with g' = 0.7, 0.8, 0.9 from the top to the bottom), while the 
dashed hnes refer to the ring approximation {g' = 0.5,0.6,0.7). The experimental data are also 
shown. 



point out that in the previous calculation the 2p-2h contributions in the A self-energy were 
evaluated by using a phenomenological parametrization of the pion-nucleus optical potential, 
also accounting for the available phase space for the 2p-2h states. The role played by 
the Landau parameters is different in the two approaches. In the present paper we have 
microscopically evaluated all the relevant diagrams which contribute at the one-boson-loop 
level; however, due to the very long computing times, we cannot implement the calculation 
for finite nuclei through the local density approximation used in 

The fairly large value of g', which we adopt in order to reproduce the A widths in the 
one-boson- loop approximation, deserves some comment. Indeed, beside the above mentioned 
phenomenology within the framework of ring approximation, which favoured values of g' be- 
tween 0.6 and 0.7, previous calculations of the (e, e') inclusive longitudinal response function 
within the one-boson-loop approximation employed a rather small value {g' = 0.35) of 
this parameter. We remind that the longitudinal response only concerns the scalar-isoscalar 
and scalar-isovector channels: hence neither the vr or p mesons nor g' enter the particle-hole 
interaction lines external to the one bosonic loop diagrams. Therefore, in ref. g' was 
employed only in conjunction with static mesons. In order to understand the influence in 
the present calculation of the dynamical (external) pion propagation, we show in figure ^ the 
results of a completely static calculation (also the "external" mesonic lines are static). The 
non-mesonic widths we obtain in this case are smaller than in figure ^ and the experimental 
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(Medium) 



2.0 



<u 

\ 



1.5 




1.0 



0.4 0.5 g' 0.6 0.7 

FIG. 8. Dependence of the non-mesonic width on the Landau parameter g', for g'^ = 0.4 and 
kp = 1.2 fm~^ in the static Umit. The various curves have the same meaning as in fig. |6[ 



band is generally reproduced by using smaller g' values: tipically g' = 0.5 in the one-loop 
approximation. This outcome somehow reconciles the parametrization of the short range 
correlations used in the evaluation of the hypernuclear decay rates with the one of ref. 
for the nuclear responses in the inclusive electron scattering. 

In table |T| we show the comparison between the one-boson-loop approximation (OBL) 
and the phenomenological model (PM) of ref. ^ at fixed kp. The calculations have been 
carried out in both cases with g' = 0.8 and g'^^ = 0.4. For technical reasons, the OBL 
calculation does not allows to precisely identify the partial ratios Fi and F2 which contribute 
to the total T^m = ^1 + ^2- In fact, we cannot separate in the imaginary parts of the 
diagrams (b)-(f) of fig. ^ the contributions coming from cuts on p-h and 2p-2h states, and 
hence the partial width (F2) stemming from the two-nucleon induced decay. The values 
listed in the table for F^^'" have been obtained from the total imaginary part of the diagrams 
^(b)-|^(f) [namely from the last two terms in the r.h.s. of eq. ( p.38|) ]. In this approximation, 
-pOBL _ p^.^^ [second, third and fourth terms in the r.h.s. of eq. (|3.38|) ]. As a matter of fact, 
one would expect that F2 increases with kp (and this is the case for the PM calculation), 
but the OBL results do not follow this statement. From table ^ and from the study 
of the 5f'-dependence of F^^^, the only conclusion we can draw on the two-body induced 
process in OBL approximation is that for LI fm^^ ^ kp ^ 1.36 fm^^ and 0.5 ^ g' ^ 1, 

r2/rfree = o.i-o.3. 
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TABLE III. Comparison between the one-boson-loop approximation (OBL) and the phe- 
nomenological model (PM) of ref. |j^ for g' = 0.8, g'j^ = 0.4. The decay rates are in unit of 
the free A width. 





kp = 


1.1 fm-i 


kp 


= 1.2 fm-i 


kp 


= 1.36 fm-i 




OBL 


PM 


OBL 


PM 


OBL 


PM 




0.82 


0.81 


1.02 


1.00 


1.36 


1.33 


r2 


0.22 


0.13 


0.26 


0.18 


0.19 


0.26 




1.04 


0.94 


1.28 


1.19 


1.55 


1.59 



V. CONCLUSIONS 

The A non-mesonic decay widths have been evaluated in nuclear matter within the 
polarization propagator method. The Feynman diagrams contributing to the A self-energy 
have been classified, via a functional approach, by the first order approximation of the 
bosonic loop expansion. In this scheme, the two-body induced decay ANN NNN is 
taken into account for the first time from a fully microscopic point of view. 

In order to compare the calculation with the available experimental data, we have em- 
ployed different "average" Fermi momenta for the following three mass regions: medium-light 
(A ^ 10); medium (A ~ 30 ^ 60); and heavy hypernuclei {A ^ 200). 

The only free parameters of our model are g' and (?^, which incorporate the NN and AA^ 
short range interactions, respectively. We have fixed these Laudau parameters in order to 
reproduce the observed decay widths, and, in the one-boson-loop approximation, the best 
choice has turned out to be g' = 0.8, g'j^ = 0.4. The agreement between the experimental 
widths and the theoretical evaluation is of the same quality as in ref. [0; interestingly, the 
same values of g' and give the best fit to the data in both approaches. However, we 
point out that the role of these Landau parameters is slightly different in the two cases. One 
difference, with respect to ref. , in the role played by g' and can be understood by com- 
paring equations ( p.8| ) and ( p. 381) : clearly, in the former (employed in ref. |I[) the short range 
correlations have a major influence in renormalizing the full polarization propagator Ul,t, 
while in the latter (used in the present calculation) the RPA renormalization only affects the 
Lindhard function Ui. Moreover, in the present calculation g' comes about together with 
static mesons (because of the discussed computational problems) in the lines "internal" to 
the bosonic loops and together with dynamical mesons in the lines "external" to the bosonic 
loops, while in all mesonic lines are dynamical. We also note that, by introducing a 
completely static description, the experimental data are reproduced with smaller g' values 
(~ 0.5). This finding is in fairly good agreement with the parametrization introduced in 
ref. for the BLE calculation of the nuclear response functions in the inclusive electron 
scattering, where g' only correlates static mesons, as explained in the previous section. 
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Finally, while the estimated rf^^ only contains one-body induced decays (in ring 
approximation), contributions to the total width stemming from this channel are also em- 
bodied in r^^'", but cannot be explicitely separated from the truly two-body induced ones. 
On the basis of reasonable extimates, the two-body stimulated decay is found to be size- 
able, (r2/rfi.ee = 0.1 -^ 0.3 for the whole hypernuclear mass spectrum studied), with values 
comparable with the results of ref. 

In conclusion, the present calculation satisfactorily reproduces the non-mesonic decay 
widths for 10 within a fully microscopic approach at fixed, average density, taking into 
account both one-body and two-body induced processes. 
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